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UNSTEADY AERODYNAMIC 
ANALYSES 


• Applications 

- Aeroelastic: blade flutter and forced vibration 

- Aeroacoustic: noise generation 

- Vibration and noise control 

- Effects of unsteadiness on performance 

• Requirements 

- Accuracy/efficiency 

* Realistic operating conditions 

* Arbitrary modes of unsteady excitation 

• Approaches 

- Numerical simulation / analytical modeling 


ASSUMPTIONS 


• Turbulence and transition can be modeled 

=3* Reynolds averaged, Navier-Stokes equations 

• High Reynolds number, “attached” flow 

=> Thin-layer Navier-Stokes equations, or 
Inviscid/ viscid interaction analyses 

• Small- amplitude unsteady excitations 

=> Nonlinear steady -I- linearized unsteady analyses 

• Re — ► oo =*> inviscid flow 

— Potential steady background flow =i> LINFLO 

— Uniform steady background flow => CLT 
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CONTRACT NAS3-25425 

NASA Program Managers: J. Gauntner, G. Stefko 

• Linearized inviscid unsteady aerodynamic analysis: LINFLO 

• Unsteady viscous layer analysis: UNSVTS 

• Steady, inviscid/ viscid interaction analysis: SFLOW-IVI 

• Coupled SFLOW-IVI /LINFLO analysis 


EXAMPLE 

CONFIGURATION 



14 


UNSTEADY EXCITATIONS 


S-oo(x,0 



• Far-field conditions (uniform mean flow) 

s(x, t) = .Re-fs-oo exp^K^ • x + w<)]}, £ < 

C(x, t) = RefC-oo expft/e-oo • x + wf)]}, ( < 

P/,too(x, t) = R e{p /iToo exp[— + tXic^ • x + wf)]}, £ > 


LINEARIZED INVISCID ANALYSES 


• Linearization 

P(x, f) = P(x) + Re{p(x) exp(t'wf)} + . . . 

► 

• Nonlinear BVP for steady background flow 

• Linear variable-coefficient problem for each Fourier 
component of first-order unsteady flow 

- Time independent 

- Surface conditions imposed at mean surfaces 

- Analytic far-field solutions for a, C, and p 

- Single extended blade-passage solution domain 

P(x -f- mG,f) = P(x) + Re{p(x) exp[t(wt + mcr)]} + . . , 

• Prescribed quantities: 

r B i S—oo C— oo» an d Pl^oo 
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LINFLO 


• Unsteady perturbation of a potential mean flow 

• Steady flow: V ■ (pV$) = 0 

• Unsteady velocity decomposition: v = V(tf» + 4>.) + v fi 

— p — —pD<f>/Dt 

— V ■ v R = 0 far upstream 
— D<t>./Dt = 0; (V4>. + y R ) ■ n = 0 on B m & W m 

• Entropy rotational velocity: X = Aer+®en -*xas(-* -oo 

s(x) — s_oo exp(iie_oo • X) 

v R (x) = [V(X • A- oo) + s-ooV$/2] x exp(*K_ 0O • X) 

• Unsteady velocity potential 

D(A~ 2 D<I>/Dt)/Dt - p -1 V ■ (pV<t>) = p -1 V • \pV<t>,] 
where <£, = F(A_ 00 , $) exp(iK_oo • X) 

• Surface conditions: 

- Blades: V<f> ■ n = /( r B ) 

- Wakes: [D<f>/Dt] = 0 and [V<£] • n = 0 

- Shocks: [pV4> + pV$] • n = /( rsh • n, V$); rsh • n = -[<A]/I$ n ] 

• Far field conditions: 

- <^/ iTCO prescribed; 4>a,*oo must be determined 

- Analytic far-field solutions for 4> = 4>i + 4 >r 
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NUMERICAL SOLUTION DOMAIN 


• Extended blade-passage region of finite extent in axial- 
flow direction 


sim 



NUMERICAL APPROXIMATION 


• Implicit, least-squares, finite-difference model 

{C4>) 0 « ( L4>) 0 — q 0 <t> o + — 4>o) 

m=t 1 

• Transonic differencing strategies 

• Cascade, local and composite mesh solutions 

• Direct solution procedure 

- Block tridiagonal system of algebraic equations for 
subsonic flow 

- Block pentadiagonal system for transonic flow with 
fitted shocks 
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AERODYNAMIC RESPONSE AT A 
BLADE SURFACE 


• Surface pressure (transonic flow): 

P(t b , t) = P(r B ) + Re{p B (r B ) exp(iwt)} + Y Ps<u( t b- <)+■•• 

n 

I 

• Blade motion: r B (x) = £(5,R,(x) 

i=l 

• Unsteady airloads: 

q, = 0 q,,rdr = -0 (P^r xe t + p B n - ^> S hJP]n] • R,dr 

• Work per cycle / pressure- displacement function 

W c =0 = 0 w(t )dr = 7r0 Im{<5*7 i r }dr = 7rlm{£ <5* 9. } 


EXAMPLE RESPONSE PREDICTIONS 

• Compressor exit guide vane (EGV): © = 15 deg, G = 0.6 

- Thick, highly-cambered NACA 0012 airfoils 

- Subsonic flow: M_ 0 0 = 0.3, H-co = 40 deg 

- Vortical excitation: u = 10, a — -2k 

- Acoustic excitation from downstream: 

UJ = 10, < 7=0 

• High speed compressor cascade: © = 45 deg, G = 1 

- Cambered NACA 0006 airfoils 

- Subsonic flow: Af-oo — 0 7, H-oo = 58 deg 

- Transonic flow: M_oo = 0.8, fl-loo = 55 deg 

- SDOF blade motions: Si — (1,0), w = 1 

• Linear/ nonlinear result comparisons 

- NGUST analysis (Navier-Stokes) 

- NPHASE analysis (Euler) 
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COMPRESSOR EXIT GUIDE VANE 

NGUST Computational Grid 



COMPRESSOR EXIT GUIDE VANE 


Steady surface pressure coefficient 



19 











VORTICITY WAVE IN AN EGV CASCADE 


FIRST HARMONIC UNSTEADY PRESSURE DIFFERENCE 
VR,-oa = (0.05g, 0), a = — 2ir, u = 10.0 



VORTICITY WAVE IN AN EGV CASCADE 


FIRST HARMONIC UNSTEADY PRESSURE DIFFERENCE 


vr-oo = (0.05g, 0), a = —2n, u> = 10.0 


— Vt-oo = 0 . 05 ? 

yjt,-co = o.io 2 

- * Vi ,-co = 0.25 ? 

VR.-00 = 0.50 2 
VISCOUS SIMULATIONS 
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COMPRESSOR EXIT GUIDE VANE 


Unsteady Pressure Response 
p+oo = ( 0 . 04 , 0 ), u> = 10 . 0 , a = 0.0 


Linearized Invisdd (LINFLO) Navier-Stokes (NGUST) 



FIRST HARMONIC UNSTEADY PRESSURE DIFFERENCE 
P/.+oo = (0.04P, 0), <7 = 0, w = 10.0 



•0.0 0.2 0.4 2cA(J 0.6 0.8 1.0 0.0 0.2 0.4 X J0.6 0.6 1.0 

REAL IMAGINARY 
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EXIT ACOUSTIC WAVE IN AN EGV CASCADE 


FIRST HARMONIC UNSTEADY PRESSURE DIFFERENCE 
P/,+oo = (0.04P, 0), <7 = 0, uj — 10.0 


■ P/, +co — 0.04P 

Pz+oo = 0.12P 

- - - - Pl % oo = 0.20 P 



NACA 0006 CASCADE 

NPHASE Computational Grid 
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HIGH SPEED COMPRESSOR CASCADE 


Surface Mach Number Distributions 
Potential, Euler 


Subsonic flow 


Transonic flow 



HIGH SPEED COMPRESSOR CASCADE 

Pressure Displacement Function Distributions for 
Torsional Blade Vibrations at a = 2 deg , w = 1 

Linearized Analysis (LINFLO) 

Nonlinear Euler Analysis (NPHASE) 

Subsonic Flow Transonic Flow 
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HIGH SPEED COMPRESSOR CASCADE 


Work per Cycle versus Interblade Phase Angle for 
Torsional Blade Vibrations at a = 2 deg, u> = 1 


□ Linearized Analysis (LINFLO) 

Nonlinear Euler Analysis (NPHASE) 

Subsonic Flow Transonic Flow 



INVISCID /VISCID INTERACTION 
ANALYSES 


• High Reynolds Number Flow 

• Inviscid region: Euler or potential flow equations 

- Surface conditions modified to account for viscous 
displacement effects 

• Viscous region: Prandtl’s equations 

- Direct solution: P —* 8 

- Inverse solution: 8 — ► P 

• Inviscid viscid interaction law 

- Weak interaction sequential solution, pressure 
determined by inviscid flow 

- Strong interaction =$» simultaneous solution, pressure 
determined by inviscid and viscous flows 
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CASCADE FLOW WITH LOCAL 
REGIONS OF STRONG INTERACTION 





Leading edge 
separation 


/ \ 

M<1 /M>1 / M<1 


Shock /boundary -layer 



Strong- 

interaction 

regions 


V-co 


interaction 


Trailing-edge 

separation 


Voo 


Weak interactions 


SFLOW-IVI: INVISCID REGION 


• Field equation 

- pV$ = 0 or ^ 2 V 2 $ = • V(V$) 2 /2 

• Surface b.c.’s account for viscous displacement effects; 
i.e., 

-Blades: V$ • n| s = p~ l d{p t u t 6) [ ds 
-Wakes: [V<t>] • n + | w = p- l d(p e u t 6 w )lds 

• Inlet flow conditions prescribed 

• Exit flow conditions determined by Kutta cond. & 
global mass conservation 
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SFLOW-IVI: VISCOUS REGION 


• Classical Viscous-Layer Eqs. (Boundary layers & Wakes) 

- Weak interaction: 

specify du e /ds — + calc. <5* (direct) 

- Strong interaction: 

specify m = p e u t 6 * — » calc. u t (inverse) 

• Turbulence and transition 

- Algebraic eddy-viscosity model 

* Blade: Cebeci-Smith w/separation modification 

* Wake: Chang, et al 

- Instantaneous transition 

• Solutions in terms of Levy Lees variables 


SFLOW-IVI: COUPLING 
PROCEDURE 
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NUMERICAL RESULTS 


• Two Cascade Configurations 

— Compressor exit guide vane (EGV) 

- High-speed compressor cascade 

• Effect of Varying Re 

• Comparison with Navier-Stokes solutions 

• Incidence Angle Study (EGV) 


COMPRESSOR EXIT GUIDE VANE 


Effect of Varying Re: Re = 10 5 

Re = 10 6 
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COMPRESSOR EXIT GUIDE VANE 


Comparison with Navier-Stokes Solution: Re = 10 6 

IVI 

N-S 
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COMPRESSOR EXIT GUIDE VANE 


Comparison with Navier-Stokes Solution: Re = 10 6 

IVI 

N-S 



0.0 0.25 0.50 0.75 1.00 

x 


29 






COMPRESSOR EXIT GUIDE VANE 


Loss Parameter, Q, & separation point location, 
x, ep , versus Inlet Flow Angle: 

Re = 10 6 , = 0.3 




COMPRESSOR EXIT GUIDE VANE 


Streamlines in Trailing-Edge Region: Re = 10 6 






HIGH SPEED COMPRESSOR CASCADE 


Comparison with Navier-Stokes Solution: Re = 10 6 

I VI 

N-S 




HIGH SPEED COMPRESSOR CASCADE 


Comparison with Navier-Stokes Solution: Re = 10 6 

IVI 

N-S 



x 
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GOAL: UNSTEADY IVI ANALYSIS 
FOR AEROELASTIC APPLICATIONS 


• High Re unsteady cascade flows 

• Inviscid region 

- Nonlinear steady (SFLOW) => $ 

- Linearized unsteady (LINFLO) =>• s,v R ,<f> 

note: V = + Re{[V(<f> + <t>') + V R ] exp(iu;f)} 

- Surface conditions 

* Blades: (V - h) • n = f B {$} 

* Wakes: [V] • n = fw{$} 

• Viscous region 

- Unsteady viscous layer analysis UNSVIS 

- UNSVIS is a direct, time marching solution proce- 
dure 

• Inviscid/ viscid coupling 

- Procedure must be developed for unsteady flows 

• Issues 

- Must modify UNSVIS to deal with moving blades 

- Matching of inviscid and viscid solutions for s and 
vr excitations 

- Need inverse unsteady viscous layer calculation 

- Inviscid/viscid coupling =$>• long computer run times; 
unless 

* 6 « 6 + 6exp(iu>t), i.e., linearization, or 

* Integral boundary layer calculation 
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INTERMEDIATE STEP: COUPLED 
SFLOW-IVI/LINFLO 

• Effects of strong steady interactions on unsteady 
pressure response 

• Assumptions 

- £(x, t ) = 6(x) + S(x, t ) 

- Strong steady inviscid/ viscid interaction 

- Weak unsteady interaction 

• SFLOW-IVI will provide steady background flow infor- 
mation for LINFLO calculation 

- Unsteady surface pressure determined by linearized 
inviscid calculation 

- Unsteady viscous layer determined by direct 
solution procedure 
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suction 



Figure 1: LINFLO results for EGV cascade undergoing torsional vibration 

(a = (1,0), cr = 0 deg, u> = 1); fl_oo = 40 deg, M_ 0 0 = 0.30: ( ) 

inviscid; ( ) viscous, Re = 10 6 . 



Figure 2: LINFLO results for EGV cascade undergoing torsional vibration 

(a = (1,0), <7 = 0 deg, u = 1); Q_oo = 54 deg, M_ oa = 0.30: ( ) 

inviscid; ( ) viscous, Re = 10 6 . 
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suction 



Figure 1: LINFLO results for HSC cascade undergoing torsional vibration 

[a = (1,0), a = 0 deg, u = 1]; = 55 deg, M_oo = 0.70: ( ) 

inviscid; ( ) viscous, Re = 10 6 . 



Figure 2: LINFLO results for HSC cascade undergoing torsional vibration 

[a = (1,0), a = ISO deg, u - 1]; — 55 deg, M- 0 0 — 0.70: ( ) 

inviscid; ( ) viscous, Re = 10 6 . 
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CONCLUDING REMARKS 


• Linearized unsteady aerodynamic analysis: LINFLO 

— Realistic 2D flow configurations 

— Arbitrary modes and frequencies of excitation 

— Efficient prediction of unsteady pressure response 

• Steady inviscid/ viscid interaction analysis: SFLOW-IVI 

— 2D cascade flows 

— Local strong inviscid/ viscid interactions 

— Efficient: CPU < 5 min 

— Robust: wide range of operating conditions 

• Future work 

— Transonic/supersonic gust response analysis 

— SFLOW-IVI /LINFLO coupling 

— Unsteady IVI analysis 



